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DETECTION METHODS AND SYSTEM USING SEQUENCED TECHNOLOGIES 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims priority of U. S. Provisional 
Application 60/425,615 filed on November 11, 2002, which is 
herein incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] This invention relates generally to the detection of a 
presence of a characteristic in an object, and, more 
particularly to the application of sequenced technologies to 
detection of a presence of a characteristic in an object. 

[0003] Detection of the presence of a characteristic in an 
object is a common function in a variety of applications. 
Exemplary applications, while not being an exhaustive listing, 
include detection of a particular disease at a location in a 
patient, detection and tracking of targets in military 
applications, detection of obstacles and objects in robotics and 
in automotive applications such as self closing sunroofs or 
doors . 

[0004] An example of detection of the presence of a 
characteristic in an object is the detection of an undesired 
substance within an object, such as detection of explosives 
inside objects. Recently there has been increased awareness of 
the potential for large-scale introduction of explosives to 
create chaos or to harm an intended set of victims. There are 
currently a variety of technologies available to detect the 
presence of explosives. However, systems utilizing a single 
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technology do not effectively measure all explosive types or 
meet the desired accuracy. In general, systems utilizing a 
single technology do not usually meet the desired conditions for 
the detection of a presence of a characteristic in an object. 
Combining technologies will not always result in improved 
performance. Data Fusion techniques can be used to provide 
guidance in the combination but such techniques require large 
computational resources and/or learning strategies. 

[0005] There is a need for simpler strategies for 
interconnecting detection systems, where each system utilizes a 
single technology, hereinafter referred to as sub-systems, and 
for combination systems designed according to that simpler 
strategy. 

Background on Explosive Detection Technologies 

[0006] Since an important example of detection of the 
presence of a characteristic in an object is the detection of 
explosives inside objects, in order to better understand the 
descriptions of systems and methods for the detection of 
explosives inside objects, a background discussion of explosive 
detection technologies is provided herein below. 

Explosive detection technologies 

[0007] The Explosive Detection Systems (EDS) considered are 
based on X-ray technologies, Neutron-Based Bulk Explosives 
Detection technologies, Other Nuclear (non-neutron) Bulk 
Explosive technologies, Nuclear Quadrupole Resonance 
technologies, and Trace Detection technologies. 
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X-Ray technologies 

[0008] X-ray Imaging, Dual Energy X-ray Imaging, Computed 
Tomography (CT), Dual Energy CT, backscatter radiography, and 
diffraction radiography all rely on x-rays penetrating into 
objects being transmitted, absorbed, scattered, or diffracted. 
Each of these processes has been exploited as an explosive 
detection system (EDS) technology. 

X-ray Imaging 

[0009] Projection x-ray imaging uses a conventional x-ray 
source exposing objects in tubs or sacks. The portion of the x- 
rays that pass through the objects are detected and converted to 
digital data and a visible image. The attenuation of x-rays is 
described by 

(1- Io/Ii = e ( ^ x) 

where Io is the output x-ray intensity, Ii is the input x-ray 
intensity, e is the natural log base, \x is the linear x-ray 
attenuation coefficient (1/cm), and x is the x-ray path length 
(cm). The [xx is actually the sum of the linear x-ray attenuation 
coefficients of each material in the x-ray path multiplied by 
the path length of each material. Each linear x-ray 
attenuation coefficient is a function of the material density 
and effective atomic number (Z). Many explosives have a density 
about 1.6 g/cm 3 and similar materials (N, 0, C, and H) so 
explosives tend to look similar in x-ray images. Unfortunately, 
many other materials have linear x-ray attenuation coefficients 
similar to explosives, including organics and plastics. A thin 
high density material can look the same as a thick low density 
material. Although projection x-ray imaging is not used for 
explosive detection, x-ray imaging can be used for identifying 
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shapes and objects related to bombs such as detonators, timers, 
batteries, and wires. Although such objects can be hidden using 
an x-ray absorbing material, such as lead, to block the x-rays, 
the lead shield would be visible in x-ray images and would cause 
a lack of penetration alarm. 

Dual Energy X-rays Imaging 

[0010] The linear x-ray attenuation coefficient of each 
material is a function of the x-ray energy, density, and 
effective atomic number. If x-ray image data is acquired with 
two or more x-ray energies (KV's) then for each picture element 
-(pixel) the density and effective atomic number can be 
calculated. Using the density and effective Z gives improved 
ability to discern explosives from similar materials. The 
density and effective Z are still averaged over the entire path 
of the x-rays, so a thin high density material can look the same 
as a thick low density material. 

Computed Tomography (CT) 

[0011] CT data is acquired by measuring the x-ray attenuation 
thru the object being examined from many different positions and 
angles. The x-rays can form a cone beam, fan beam, or pencil 
beam. The x-rays can be detected as an area, a strip, a line, 
or a single point. An x-ray source emits x-rays that are 
attenuated as they pass through the object being examined. The 
intensity of the x-rays exiting the object are measured by the 
detectors, digitized, and fed to a computer along with position 
information. The computer uses the x-ray three dimensional 
vector information to determine the location of each voxel 
(volume element) that each x-ray passed through. The process 
called back projection adds part of each x-ray measurement to 
each of the voxels (volume elements) that the x-rays passed 
through. The result is a 3 dimensional (3D) reconstruction of 



the object being examined. The linear x-ray attenuation 
coefficient , \k, is estimated for every voxel. More x-rays and 
more angles passing thru a voxel give better estimates of u. CT 
provides a 3D array of measured x-ray attenuation coefficients 
as a function of position.. As with projection x-ray, the x-ray 
attenuation coefficient is directly related to the material in 
the voxels. CT has an advantage over projection x-ray in that 
instead of averaging over the entire x-ray path, there is an x- 
ray attenuation coefficient for each voxel. Smaller voxels make 
it less likely that more than one material is in the voxel and 
likelihood of correctly interpreting the material is improved. 
The CT x-ray source and detector can be mounted on a gantry, 
which rotates around the object being examined. Similarly, the 
CT x-ray source and detector can be fixed where the object being 
examined is rotated. Typically, a fan beam is used with a line 
detector on a gantry to gather slice data and the object is 
translated on a conveyor to acquire multiple slices. It should 
be noted that the CT X-ray source could also be a cone beam 
source. Furthermore, it should be realized that the above 
embodiments are sample embodiments and that other embodiments 
are possible. If the conveyor is fixed during acquisition the 
result is a slice and if the conveyor is moving the result is a 
helical scan. 

[0012] Alternatively, the CT can use multiple x-ray sources 
and multiple detector arrays to eliminate the need for a 
rotation. Fix source CT has advantages of less moving parts, 
higher speed, higher mechanical reliability and disadvantages of 
lower electrical reliability, more electrical parts, less 
accurate reconstructions due to fewer measurements. Each 
configuration requires specialized algorithms for back 
projection and reconstruction. All CT configurations measure 
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the x-ray attenuation coefficient for each voxel, however some 
systems produce small voxels covering the entire volume while 
some have space in between the voxels and some systems only look 
at the parts of the volumes that were found to be suspect by 
other technologies. All are types of CT, but their capabilities 
vary widely. 

Dual Energy CT 

[0013] Two sets of CT data acquired using two different x-ray 
energies can be used to generate reconstructions for the two 
energy levels. The linear x-ray attenuation coefficients for 
the two energies for each voxel can be used to solve for the 
density and effective atomic number of each voxel. The density 
and effective atomic number can give a more accurate indication 
of the material in the voxel than just x-ray attenuation 
coefficients of single energy CT. Dual Energy CT (DECT) is more 
effective than dual energy projection x-ray because multiple 
materials are less likely to be found in a voxel than in an 
entire x-ray path of projection x-ray images. Multiple energies 
can be used to further enhance the data quality, but not by 
much. X-ray sources generally emit wide spectrums of x-ray 
energy that tend to overlap and blur the value of more than two 
energies. The compromise with DECT is that widely separated 
energies work better, but widely separated means that one energy 
must be lower energy and the lower energy may not penetrate the 
object being examined as well. The number of x-ray photons may 
be too low to give a useful improvement in the probability of 
detection and probability of false alarm for explosive detection 
in some objects. 

Backscatter Radiography 

[0014] When x-rays enter a material some of the x-rays are 
absorbed by materials and reemitted (scattered) at lower energy 



in all directions. Those x-rays returning to near the x-ray 
source are called backscatter. A collimated pencil beam x-ray 
source is used to determine the location doing the 
backscattering. Scanning the x-ray pencil beam can produce an 
image of the backscattered x-rays. Backscatter images are 
particularly good at imaging lower density materials such as 
organics, plastics, and explosives. Backscatter imaging can not 
distinguish between the material types. Scanning pencil beam 
imaging systems, such as those used in backscattering imaging 
systems, may simultaneously produce a transmission or projection 
x-ray image which is useful in for identifying shapes and 
objects related to bombs such as detonators, timers, batteries, 
and wires. Backscattering is a self -filtering process so the 
farther into the object being examined the x-rays travel, the 
less sensitive backscatter image is. 

Diffraction Radiography 

[0015] When x-rays entering a material encounter crystalline 
materials, the x-rays are bent or diffracted in a direction 
dependent on the material. Some explosives components are 
crystals and can be detected by the diffraction angle of the x- 
rays. X-rays are diffracted by many other common materials 
including metals, sand, gems, medicine, some papers, and photos. 
Thick materials give wide diffuse diffraction angles making it 
hard to distinguish one material from another. Simply detecting 
crystalline materials gives insufficient explosive detection 
specificity. This is another technology that is useful in the 
lab and with special test objects, but may be difficult to use 
in a real application. 
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Neutron-Based Bulk Explosives Detection Technologies 

[0016] The basic principles behind the neutron-based 
Explosives Detection Systems (EDS) are the following: (1) a beam 
of neutrons enters the package to be inspected, (2) the neutrons 
interact with nuclei present in the explosives, (3) as a result 
of the interaction, either (a) new radiation, characteristic of 
the target nucleus, is produced, and/or (b) a deflection or 
absorption of the neutron beam occurs, (4) the new radiation 
and/or the neutron beam exiting the target is detected, and (5) 
computer analysis of the outgoing particles is used to deduce 
the presence of the explosive. All of the neutron-based 
techniques try to detect explosives by detecting and quantifying 
the atomic components of explosives, nitrogen (N) , oxygen (0), 
carbon (C) and hydrogen (H) . The main signature for explosives 
is the relatively high concentration of nitrogen and secondarily 
oxygen, and the relatively low concentration of carbon and 
hydrogen, as compared to common materials. The various neutron- 
based technologies differ in the radiation (referred to in step 
3, above) they use to deduce the presence of one or more of N, 
0, C and H, the properties of the neutron beam whose interaction 
produces the radiation, and the method, if any, used to localize 
(i.e., image) those materials within the package being 
inspected. In the following, a brief description of the 
theoretical basis of each of the neutron-based technologies is 
provided, starting with the most mature technology and ending 
with the least. 

Thermal Neutron Analysis (TNA) 

[0017] In TNA the detected radiation is a gamma ray. The 
neutron interaction producing that gamma ray is the capture of 
the neutron by the nucleus; for example, the isotope 14 N captures 
a neutron to become 15 N. Although 15 N is a stable isotope, it is 
initially produced in a very highly excited state, which 
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subsequently (within 10" 17 s of the capture) decays to its stable 
state, producing a number of gamma rays, the most energetic of 
which has energy of 10.8 MeV. The probability of neutron capture 
is (usually) highest for the slowest neutrons; the slowest that 
a neutron inside a target can be is when its energy is the same 
as the average thermal energy of the target, hence the name TNA. 
The cross section for thermal capture on 14 N is 75 mb 1 , on 16 0 0. 19 
mb, on 12 C 3.53 mb, and on X H 333 mb. (The cross section is a 
measure of the likelihood that a process take place; the above 
numbers indicate that for an equal number of atoms, a thermal 
neutron is 75 mb / 0.19 mb = 394 times more likely to be 
captured by a 14 N nucleus than by an 16 0 nucleus, and 333 mb / 75 
mb = 4.4 times as likely to be captured by a l E nucleus than by a 
14 N nucleus.) These cross sections indicate that TNA is primarily 
sensitive to the N component of explosives and not to the 0 
content. Thus the presence of relatively small quantities of N- 
rich substances such as Polyurethane or Melamine, or large 
quantities of substances containing small-to-moderate 
concentrations of N, such as Orion, silk and wool, can trigger a 
false alarm in a TNA system 

[0018] The neutron source for TNA could be a radioactive 
source or an accelerator based source (the term Electronic 
Neutron Generator, or ENG, is sometimes used in the industry to 
refer to accelerator-based neutron sources). The most suitable 
radioactive source for TNA is the fission source Californium 
252 , 252 Cf, which produces 2.30xl0 6 n/s per microgram. The energy 
distribution of the emitted neutrons peaks between 0.5 and 1.0 
MeV and extends to as high as 10 MeV, with the intensity 
dropping almost exponentially with increasing neutron energy. 



1 For N, thermal capture produces a 10.8 MeV gamma ray only 14% of the time, so the effective cross section for 
that signature is 0.14*75=1 1 mb. 
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Thus the size of the source is determined by the encapsulation 
requirements and the amount of moderator needed to bring down 
the energy to thermal levels. The half -life of 252 Cf is 2.6 years. 

[0019] Accelerator based neutron sources for TNA make use of 
either the D-D (deuteron-deuteron) reaction: 2 H + 2 H 3 He + n, 
or the D-T (deuterium-tritium) reaction: 2 H + 3 H -» 4 He + n. The 
deuteron (nucleus of the deuterium heavy isotope of hydrogen) is 
accelerated by a potential of 100-300 kV and the accelerated 
beam usually impacts a deuterated or tritiated target, such as 
titanium deuteride, or titanium triteride. The D-D reaction 
releases 3.26 MeV of energy and the D-T reaction releases 17.6 
MeV of energy. Because the released energy is much larger than 
the beam energy, all outgoing neutrons leave with about the same 
energy, approximately 2.7 MeV for the D-D reaction, 14 MeV for 
the D-T reaction. A 1 mA beam of deuterons will produce about 10 9 
n/s for D-D reaction and 10 u n/s for the D-T reaction; somewhat 
smaller yields are produced from compact sealed-tube neutron 
generators. Because of the depletion of the target, sealed tube 
generators have operating life times of less than 2000 hours. As 
in the case of a radioactive source, the neutrons are produced 
isotropically and hence one can make use of only those neutrons 
that end up in the direction of the target to be inspected. 

[0020] The gamma rays produced in TNA are usually detected 
with Nal(Tl) or BGO crystals. The gamma-ray spectrum usually 
consists of a full-energy photo peak plus two escape peaks plus 
a continuum Compton background. For the energies encountered in 
TNA and FNA (see below) typical sizes are 4"x4"x4" per crystal; 
while larger size crystals do exist, larger sizes result in 
higher count rates and in TNA could easily reach a MHz, which is 
probably the upper limit for single pulse counting with Nal . One 
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then employs as many crystals around the package being inspected 
as geometry and cost allow. 

[0021] Other scintillating crystals, such as BGO, have a 
higher efficiency than Nal, but have a lower energy resolution 
and a higher cost. Given the potential complexity of the gamma- 
ray spectra, the better energy resolution of Nal crystals favors 
their use. 

Fast Neutron Analysis (FNA) 

[0022] As in TNA, the detected radiation signaling the 
presence of an explosive is a gamma ray characteristic of one or 
more of its atomic components. The neutron interaction producing 
the gamma ray is, however, not thermal capture, but inelastic 
neutron scattering. In this process an energetic (i.e., fast, 
hence the name FNA) neutron collides with a nucleus and 
transfers some of its energy to that nucleus. The nucleus then 
decays back to lower excited states giving off one or more 
characteristic gamma rays . The minimum neutron energy needed to 
excite the first gamma-decaying state of 12 C, 14 N, and 16 0 is 4.4 
MeV, 2.3 MeV, and 6.2 MeV, respectively. One of the advantages 
of FNA over TNA is its sensitivity to 16 0 and 12 C, in addition to 
its sensitivity to 14 N. This sensitivity comes about because the 
cross section for inelastic scattering from those atoms is much 
larger than the cross section for thermal neutron capture. For 
14.8 MeV neutrons (i.e., those produced by the D-T reaction) the 
cross section to excite the above states is 200 mb for 12 C, 70 mb 
for 14 N, and 100 mb for 16 0. However, since the minimum neutron 
energy needed to excite the oxygen is 6.2 MeV, one cannot use 
the radioactive source 252 Cf as the neutron source; one must use 
either the D-T reaction, which produces 14.8 MeV neutrons for an 
accelerator that provides the minimum of -200 keV of energy to 



the accelerated deuteron, or the D-D reaction, with an 
accelerator that provides at least 5 MeV to the deuteron. (Other 
means of producing fast neutrons exist, but only the methods 
employed in machines produced purely for explosives detection 
using neutron-based technology are described herein.) 

[0023] The gamma-ray detectors used in FNA are of the same 
type used in TNA, (mainly Nal) since the energies of interest 
are of the same order, with FNA perhaps affording the use of 
slightly smaller detectors since the maximum energy of interest 
(6.13 MeV from 16 0) is smaller than that encountered in TNA (10.8 
MeV from capture on 16 N) . 

[0024] The improvement in FNA over TNA comes about because of 
the enhanced sensitivity to 16 0, whose density is highly 
correlated with that of N in explosives (both 0 and N content is 
high in explosives). However, it is possible to get a false 
identification if it happens that there are two (or more) 
objects in the package being inspected one with high N content 
and the other with high 0 content, with the ratio being close to 
that of an explosive. To reduce the probability of such false 
positives, it is desirable to add imaging capability to the 
above techniques, i.e., is it is desirable to add the ability to 
localize the position of the N and 0 signatures within the 
package. The theoretical basis for variations of the FNA and TNA 
techniques, designed to add imaging capability and/or improve 
their sensitivity, is presented herein below. 
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Short Pulse FNA (SPFNA, also known as nanosecond pulse FNA or 
simply as PFNA) . 

[0025] This technique is based on FNA. The deuteron beam 
producing the neutrons through its interaction with a deuterium 
target is pulsed, i.e., the beam particles are bunched together 
(by electromagnetic means) into pulses each about 1 ns (10~ 9 s) 
in duration and the pulses repeated about 1-5 MHZ, i.e., the 
time between one pulse and the next is about 200 — 1000 ns . The 
reason for this short pulsing is to determine where the neutron 
interacted in the package being inspected, and hence to localize 
the source of the gamma rays (and hence the source of the 0 and 
N). The localization along the direction of travel of the 
neutron is deduced as follows: the electronics that produce the 
deuteron pulse in the beam generate a fast time signal with each 
pulse. This time signal starts an electronic timer; when a gamma 
ray is detected in the Nal (say) the generated pulse is used to 
read the timer. (In practice the timer is started on the Nal 
signal and stopped on a delayed version of the beam pulse 
signal.) The measured time is linearly proportional to the 
distance traveled by the neutron and hence gives the position 
along the neutron travel direction where the scatterer is 
located. (The effect of the additional time taken by the gamma 
ray to travel from the neutron interaction site to the detector 
is discussed below. ) The position in the two directions 
perpendicular to the neutron travel direction is determined by 
sufficient collimation of the neutron beam and by rasterizing 
the beam along one of the two directions (say the vertical 
direction) while a conveyor belt carrying the package being 
inspected and moving at a known speed rasterizes the object 
along the horizontal direction. 
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[0026] The position resolution (along the beam direction) is 
determined by the speed of the neutron, the duration of the 
pulse, and the time resolution of the Nal-electronics 
combination. For an 8.4-MeV neutron, the speed is 2.8 cm/ns; for 
a Nal detector the time resolution is at best 1 ns; when 
combined in a statistical way with the 1 ns pulse width the 
resulting time resolution is probably at best 1.4 ns, leading to 
a position resolution of about 2.8 cm/ns * 1.4 ns = 4.0 cm. (= 
1.6"). For a 14.8 MeV neutron, such as would be produced by a 
low- voltage D-T accelerator, the neutron speed is 3.7 cm/ns and 
the corresponding resolution would be 2.1". The resolution in 
the direction perpendicular to the beam direction is determined 
primarily by the collimation of the neutron beam and other 
geometrical considerations. Because of the much larger speed of 
gamma-rays (the speed of light = 30 cm/ns), the additional time 
taken by the gamma ray to travel from the point of interaction 
to the detector is small compared to that taken by the neutron 
to travel within the package; nevertheless, with a knowledge of 
the position of the detector, its contribution can be taken into 
account with sufficient accuracy so as not to degrade the 
position resolution. 

[0027] One factor contributing to the degradation in position 
resolution is scattering of the neutron prior to its gamma- 
registering interaction; the probability of that interaction 
depends on the size and content of the package being inspected. 

Long Pulse FNA (LPFNA, microsecond pulse FNA or Gated FNA, or 
gated FNA/TNA) 

[0028] Like in SPFNA, the deuteron beam producing the 
neutrons is pulsed, however the pulse duration here is on the 
order of several microseconds (typically 10-50 \is) , and the 
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separation between pulses is on the order of 100-1000 \is. The 
reason for the pulsing here is not to obtain position 
information (the pulse duration and the resulting spread in 
interaction positions is far too large to provide useful 
imaging), but rather to improve the sensitivity and accuracy of 
the gamma-ray measurement. This is done by separating those 
gamma counts that arrive during the pulse (i.e., with the beam 
on), from those that arrive in between pulses (beam off). The 
former are then largely due to fast neutron interactions (i.e., 
close to what would be observed in an idealized FNA 
measurement), and the latter largely due to thermal neutron 
interactions (i.e., close to what would be observed in an 
idealized TNA measurement). In principle one can also employ a 
"pure" gated TNA technique (GTNA), in which the neutrons are 
thermalized with a moderator before they get to the package and 
then only count gammas during the beam off period. It should be 
clear that in either case, the pulsing is still too short to be 
done by mechanical means and hence only electronically generated 
neutrons (i.e., accelerator-based) can be used. 

Neutron Resonance Attenuation (NRA) 

[0029] In this technique a ns-pulsed beam of neutrons, with a 
broad energy spectrum (0.5-5 MeV) enters the package to be 
inspected from one side, and the neutrons that pass through the 
package are detected on the other side. The elements making up 
explosives (N, 0, C, H) have large scattering cross sections at 
specific energies (that is resonances in the cross section) and 
hence , if present, would scatter away neutrons of that energy 
from the beam (i.e, the beam component at that energy would be 
attenuated) . Therefore by looking for depletions in the 
transmitted neutron beam at characteristic energies the presence 
of explosives constituents could be deduced; the depth of the 
depletion provides a quantitative measure of the total areal 
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density (density x thickness) of the scattering element along 
the beam direction. By rasterizing the beam (and/or the target), 
a two-dimensional shadowgram of the elements of interest can be 
formed. The neutron energy is usually measured with a Time-of- 
Flight (ToF) technique: the beam pulse starts a timer and the 
neutron detector pulse stops the timer. The measured ToF then 
gives information on the velocity of the neutron (and hence its 
energy), since the distance between the neutron source and the 
detector is fixed. The sensitivity of the technique depends, 
among other factors, on the energy resolution, which, in turn, 
depends on the length of the flight path, with longer flight 
paths giving higher resolution. Typical flight paths are 2—10 
m long. 

Associated Particle Imaging (API) 

[0030] This technique uses the neutrons to interrogate the 
target in essentially the same way FNA does: a fast neutron 
excites the N, 0 and C nuclei and a subsequent gamma 
characteristic of these nuclei is detected. The imaging is done 
in a different way than in SPFNA, though. The reaction used to 
produce the neutron is the D-T reaction: 2 H + 3 H -» 4 He + n. In 
this reaction the direction of the alpha particle (the nucleus 
of 4 H) is correlated with the direction of the neutron (for the 
low beam energies typically used in the D-T reaction the alpha 
particle and the neutron are emitted essentially back-to-back). 
Therefore if the alpha particle is detected at a known direction 
(say by using a position sensitive detector and/or an array of 
detectors each with a small angular size) then the direction of 
the associated neutron is known. The position of the N, C, 0 
along the (now known) direction of the neutron is obtained from 
the time difference between the alpha particle and the gamma 
ray. 
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[0031] While API would allow the use of a simpler accelerator 
than used in SPFNA, and would also allow a larger area of the 
target to be inspected at a time, the high count rates in the 
alpha detectors necessitate the use of lower beam currents and 
hence longer inspection times. While there may be ways to 
ameliorate the alpha count rate problem, for example by using a 
large number of small detectors each with a lower count rate, 
such systems have not, to our knowledge, been developed 
commercially yet. 

Other Nuclear (Non-Neutron) -Based Bulk Explosives Detection 
Technologies 

[0032] An overview of other nuclear-based techniques for 
explosives detection that do not rely on the generation and 
interaction of neutrons follows. The theoretical basis for two 
technologies that are under active commercial development, Gamma 
Resonance Technology (GRT) and Gamma Radiography, is presented 
herein below. 

Gamma Resonance Technology (GRT) 

[0033] As in Fast Neutron Analysis ( FNA ) , GRT detects 
explosives by detecting the presence of one of the principal 
components of explosives, nitrogen. The nitrogen is detected by 
probing the package with a gamma ray whose energy is uniquely 
tuned to have an enhanced interaction probability with nitrogen. 
Further, the gamma-ray beam is produced in such a way that only 
nitrogen along a specified path interacts with the beam, thus 
allowing imaging of the nitrogen. The way these unique gamma 
rays are produced is as follows: A beam of protons with an 
energy of 1.75 MeV impinges on a 13 C target. Some of the protons 
are captured by the 13 C nuclei to form 14 N in an excited state at 
9.172 MeV. About 5% of the time that excited state decays by 
emitting a gamma ray with an energy of 9.17 MeV (the remaining 
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95% of the time the state decays back into a proton plus a 13 C 
nucleus). Now the laws of physics, at the microscopic level, 
seem to be time reversible, so that if the process in which a 
stationary excited 14 N nucleus decays to the ground state by 
emitting a 9.17 MeV ray happens, then its time inverse, namely 
the process in which a 14 N nucleus in its ground state absorbs a 
9.17 MeV gamma ray to become an excited 14 N at 9.172 MeV 
excitation should also happen, with essentially equal 
probability. So the presence of 14 N in the package being 
inspected can be inferred by monitoring the number of 9.17 MeV 
gamma rays passing through the package. A decrease in the number 
of transmitted gammas is indicative of the presence of 14 N, 
potentially signaling the presence of an explosive. Now there, is 
a slight complication in the physics which has both negative and 
positive consequences for explosives detection. The excited 14 N 
produced by the beam is not stationary, but is moving in the 
direction of the beam (it has to carry the momentum brought in 
by the proton). Because of this motion, the gamma rays emitted 
in the forward direction (i.e., in the same direction as the 
proton beam) have more energy than can be absorbed by a nitrogen 
at "rest", while gamma rays emitted in the backward direction 
have less energy and cannot excite the 14 N to the excited state 
at 9.172 MeV. Consequently, there is a specific angle (relative 
to the proton beam direction) at which the gamma rays have just 
the right energy to be absorbed by the 14 N. This angle is about 
80. 5°. 2 The natural width of the 9.172 MeV level in 14 N is 122 eV, 
so that gamma rays with about that much spread in energy can 
still be resonantly absorbed by the 14 N; the spread in angle 
allowing that spread in energy is about 0.2°. 



2 At that angle the energy of the gamma ray emitted by the moving 14 N is 9.175 MeV, about 3 keV higher than the 
excitation energy of the resonant state in 14 N. This extra 3 keV is needed because upon absorption of the gamma ray 
the ,4 N recoils back and hence takes up 3 keV of the energy of the gamma ray. 
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The negative consequence to the fact that only gamma rays 
emitted within less than half a degree from a cone with half 
angle of 80.5° from the beam direction is that only a small 
fraction of the emitted gamma rays (about 0.4%) fall within that 
range. The positive consequence is that the narrow resonant fan 
beam of gammas can be used to image the nitrogen content of a 
package as it passes through the slice , without the need for 
very precise collimation of the gamma beam. 

[0034] There are two ways the presence of 14 N in the package 
under inspection can be deduced: 1) a decrease in the number of 
counts of 9.172-MeV gammas along the direction of the cone. To 
be sensitive to that decrease, a detector capable of 
discriminating between the resonant gammas and the non-resonant 
gammas is needed. The principle used is essentially the same one 
used in detecting the 14 N: a detector with 14 N content is used; 
some fraction of the transmitted resonant gammas will be 
resonantly absorbed by the 14 N in the detector; the 1.75-MeV 
protons from the decay of the excited 14 N to p+ 13 C are detected. 

[0035] 2) Detection of 9.1 MeV gammas resonantly scattered by 
the 14 N away from the cone. Such a detector would be placed at 
relatively backward angles from the resonant gamma fan. The 
detector need only have moderate energy resolution for 9 MeV 
gammas, since the probability of non-resonant scattering is very 
small . 

Gamma Radiography 

[0036] Gamma Radiography employs exactly the same principle 
as x ray radiography: The variation in the absorption of a beam 
of photons as an object is scanned by the beam is used to from a 
shadowgram of the object. The shadowgram is produced because 
denser objects are more absorbing than less dense objects. In 



gamma radiography the gamma rays are produced by the decay of a 
radioactive isotope, instead of an x ray tube, or an electron 
accelerator. Typical sources used are 137 Cs, which produces a 662- 
keV gamma ray, and 60 Co, which produces two gamma rays, one at 
1172 keV and the other 1332 keV. The higher energy of these 
gamma ray sources, compared to typical x ray sources, enables 
the examination of thicker targets; by the same token this 
higher penetration means that the contrast in the shadowgrams is 
relatively low. 

Nuclear Quadrupole Resonance (NQR or QR) 

[0037] Quadrupole resonance (QR, or more formally, nuclear 
quadrupole resonance, NQR) offers a highly specific method of 
detecting and identifying crystalline materials, including many 
explosives. Specially tuned pulses of radio waves in the 0.5-5 
MHz range are applied to the sample. The radio waves are 
absorbed by nuclei of atoms in crystals within the sample and 
re-emitted a short time later. These weak re-emitted radio 
waves are detected by a sensor coil in the NQR receiver. The 
radio spectrum corresponds to particular atoms in a particular 
molecular structure, so each peak in the spectrum indicates a 
unique crystalline material. 

[0038] Not all atomic nuclei respond to the QR radio pulses. 
The resonance effect is specific to nuclei with magnetic spin 
quantum numbers greater than or equal to 1, which includes 
Nitrogen-14 and other nuclei typically found in explosives. 
Another requirement for QR is that the nuclei must reside in an 
environment with a static or non-temporally averaged electric 
field gradient (EFG). What this means practically is that the 
nuclei must exist in an environment where the arrangement of 
surrounding charges remains fixed for at least a few 
microseconds. Nuclei in solid, amorphous materials could posses 



a QR spectrum, but that spectrum would tend to be broad and 
indistinct reflecting the presence of many different EFG 
environments. Nuclei in crystalline environments where atomic 
arrangements regularly repeat and the atoms don't move much will 
have narrower QR spectra as long the charge arrangement 
surrounding each nucleus gives a non-zero EFG, This describes 
the situation for many explosive and drug compounds. Nuclei in 
liquids do not possess a QR spectrum because the EFG is 
constantly changing and averages to zero much faster than 1 
microsecond. 

[0039] Although somewhat related to Magnetic Resonance 
Imaging used in medicine, QR is distinctly different in several 
ways. The most notable difference is that QR does not use a 
magnet to polarize the nuclei. QR polarization results from 
different energy states for the nuclei in the EFG. Not having 
to use a magnet is a major advantage of QR over competing 
techniques based on other magnetic resonance effects. Some of 
these advantages are lower cost, simpler design, and easier 
portability. The other difference between QR and MR I is that QR 
normally does not yield an image. A positive QR signal 
characteristic of an explosive indicates the presence of an 
explosive somewhere within the sample. Some degree of 
localization may be possible by taking advantage of the spatial 
distribution of RF field produced by specific antenna designs. 
Like MRI, QR uses RF radiation that is non-ionizing, so that it 
is generally safe for use on or near personnel. 

[0040] The RF pulses are feed into a coil surrounding or 
adjacent to the sample, thereby producing RF fields within the 
sample. Likewise, the signal pickup coil must surround or be 
near the sample. In some cases, the transmit and receive coils 
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are one in the same. Design of these coils as effective 
antennas for each NQR application and the associated electronics 
to handle and analyze the signals are major aspects of NQR 
technology development. 

[0041] Because this technology uses radio frequency fields , 
RF shielding of the item under inspection becomes a possible 
means for defeat. Separate means can be devised to detect when 
an item being inspected contains a shielded region. 

Trace Detection Technologies 

(a) Trace Detection Technologies Summary 

[0042] Intrusion into the package is required to limit false 
negative objects, i.e. not seeing a well packaged explosive, to 
an acceptable level. Once intrusion into the package is 
acceptable Trace Detection Technologies are appropriate to 
identify specific chemical explosive and appropriate emergency 
response. The various trace detection technologies are 
described in the following sections. 

(b) Trace Detection Technology Description: 

[0043] Trace Detection Technologies require physical 
transport of the sample (explosive) to the detector. Detectors 
only respond to (detect) gas phase chemical species. 
Particulate transported to a detector must go through a 
vaporization phase change prior to detection. Detection limits 
are best defined in absolute mass at the detector. Sampling, or 
transporting the chemical to the detector, is the most important 
limitation of trace technologies. 
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[0044] To illustrate the significance of packaging on 
detection explosive systems, two cases are described herein 
below. 

[0045] The technology is adequately described by breaking the 
technology into three parts; sampling, sample preparation, and 
detection. The three independent parts are detailed in Table 1. 
First, define the analytical situation or problem, and then 
select one technique from each column appropriate to maximize 
and tailor the Trace Detection Technology to the problem. 



Table 1: Trace Detection Technologies have three independent parts that are mixed and matched to achieve 

analytical goals 



Sampling Method 


Sample Preparation 


Detectors 


Air, human plume 


None or direct injection 


Ion Mobility Spectrometer 


Air, portal 


Chromatography 


Mass Spectrometer 


Air w/absorbent 
preconcentration 


Cryo-fractionating 


Infrared spectrometer 


Wipe sample 


Thermal desorption 


Electron Capture Detector 
(ECD) 



[0046] Example detector sensitivities are listed in Table 2. 
These are estimates based on ideal laboratory tests. Detector 
sensitivity is often compound (chemical functionality) specific. 
In real world applications, significantly worse detection limits 
should be expected. Much of the sensitivity lost by taking 
instrumentation out of the laboratory is regained through a pre- 
concentration step as part of the sampling. Little difference 
in sensitivity is seen as a discriminator between manufacturers 
at these levels of detection once installed in the field. 

[0047] Mass spectrometers (MS) and infrared (IR) 
spectrometers are generally considered separately from other 
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detectors as they provided chemical structural information. 
This additional dimension to the detector provides significant 
advantages over nonspecific detectors. Using MS and IR 
spectrometers, the explosive in the object can be identified. 



Table 2: Detector Sensitivity Limits 



Detector Type 


Lower 
Detection limit 

(gm) 


Linear Range 


Use 


Nonspecific Detector Technologies 


Thermal 
Conductivity 


10 A -6 


10M 


General 


Flame Ionization 


10M0 


10 A 6 


Organic 


Electron Capture 


10M2 


10 A 2 


Variable 


Chemiluminescence 


10M0 


10 A 5 


Oxides of nitrogen 


Surface Acoustic 
Wave 


10M2 


10 A 5 


General 


Cross-section 
Ionization 


20 A -6 


10 A 5 


General 


Phosphorus 


10 A -11 


10M 


Organo-P 
Compounds 


A little better than Nonspecific 


Ion Mobility 
Spectrometry 


10M2 


10M 


High Proton 
Affinity 


Specific Detector Technologies 


Mass Spectrometer 


10 A -9 


10M 


General 


Infrared spectrometer 


10 A -9 


10M 


General 



[0048] Ion Mobility Spectrometer (IMS) detector are one step 
up from a nonspecific detector. What is compromised in 
resolution is gained in speed, usability, weight and cost. 
Within limits IMS can be used to identify chemical species 
through the use of computational analysis. 
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[0049] Chromatography cannot be used as a stand-alone 
detector technology. It can be used with most detectors and 
functionally prepares the sample, separates chemical species to 
facilitate the use of nonspecific detectors, or simplify 
computational analysis of MS and IR detectors when complex 
mixtures are sampled. 

[0050] These technologies, and specifically MS/RGA with a 
thermal sample inlet, are appropriate to use when the package 
can be intruded upon and the specific chemical explosive needs 
to be identified to guide an emergency response. 

Mass Spectrometry ; 

[0051] Mass spectrometer detectors fall into two types; bench 
units primarily used in laboratory settings and require highly 
skilled operators, and mobile units. Mobile units designed for 
gas analysis are often called RGAs (residual gas analyzers), and 
are often equipped with a thermal sample inlet. The sample 
inlet serves two (2) functions. First, sample is drawn into the 
sample inlet cold to concentrate chemical species and thus 
increase sensitivity. Second, the inlet is quickly heated 
vaporizing collected chemical species for introduction to the 
mass spectrometer detector allowing measurement of a variety of 
low vapor pressure chemical species at higher sensitivity. 

[0052] The mass spectrometer detector has an electron 
ionization (about 70eV) source. As chemical species enter the 
ionization beam positive ions form. The mass to charge ration 
of these ions and the uni-molecular decomposition, or 
fragmentation, products are measured. Every molecule has a 
unique mass fragmentation pattern that can be used to separate 
and identify it from any other compound. Greater computing 
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power is needed when either more components are in the sample 
mixture or the more unknown chemical species are in the sample. 
To lessen the required computing power chemical species are 
separated or fractioned through the use of a short 
chromatographic column placed between the thermal inlet and the 
detector . 

[0053] Mass spectrometer detectors greatest strength lies in 
its ability to provide structural proof of the chemical species. 
When screening level gets down to the point where chemical 
speciation is required (specific explosive or drug), this is the 
technology to use. 

Ion Mobility Spectrometer (IMS) 

[0054] IMS is well developed and in various forms has been 
around for years. Samplings offered by the listed IMS venders 
are air flow (human plumb or portals) or wipe, sample 
preparation is direct injection and detection is by IMS using 
atmospheric pressure chemical ionization (APCI). 

[0055] IMS sensitivity is dependent upon the chemical 
species, generally proportional to higher proton affinity, and 
can range over several orders of magnitude (1CT 12 gm to 1(T 6 gm) . 

[0056] Ion mobility spectrometers use a radioactive source 
(about 15 milli-curies Ni63, some use Am241) to ionize gas 
molecules at atmospheric pressure (APCI atmospheric pressure 
chemical ionization) . Ion clusters form and drift in a field, 
drift tube, being detected as they collide with spaced detector 
elements along the tube. The time of drift (mobility) is 
roughly a measure of the cluster mass. Often referred to as a 
time-of -flight measurement. Its limited chemical speciation is 
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achieved based on the drift time or mobility in a positive 
and/or negative field. This is very oversimplified. The point 
is that IMS provides rudimentary identification of chemical 
species giving it advantage over a parallel technology, electron 
captor (ECD). It gives significantly less identification 
capability when compared to mass spectrometry (MS). Most drugs 
and explosives have high proton affinity supporting IMS high 
detection sensitivities to these classes of compounds. 

Electron Capture Detector 

[0057] Electron capture detector (ECD) is a nonspecific 
detector. ECD use a radioactive beta source to ionize a gas 
inside a cell. The current at the collector electrode is 
measured. As chemical species inter the cell they react with 
the ionized gas decreasing current flow to the detector. As ECD 
is very sensitive to electronegative chemical compounds and very 
insensitive to neutral compounds it is considered a selective 
nonspecific detector. 

Chemiluminescence 

[0058] The Chemiluninescence detector, although specific for 
oxides of nitrogen, is considered a general nonspecific detector 
with high sensitivity. 

[0059] A detector technology generally used to measure the 
concentration of oxides of nitrogen (NOx). "Nitrogen oxide (NO) 
and ozone (03) react to produce a characteristic luminescence 
with an intensity linearly proportional to the nitrogen oxide 
concentration", Model 42c Instruction Manual, Thermo 
Environmental Instruments. Nitrogen dioxide is first reduced by 
passing over a hot catalyst followed by mixing with ozone. Some 
additional specificity is achieved by adjusting the catalyst and 
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the catalyst temperature. Chemical species with high nitrogen 
content will react. Some explosives , energetic materials, and 
drugs react similarly when reacted with ozone either before or 
after the reduction catalyst. 

Olfaction 

[0060] Technologies offered under the Olfaction name are more 
a way of handling data than a specific detector. An array of 
very small simple materials that differ in chemical properties 
(i.e. solubility or polarity) will interact or react with a 
specific chemical species (TNT or cocaine or anthrax) 
differentially exhibiting a unique response pattern. The more 
elements in the array the more unique the pattern. Thirty-two 
elements are common. Chemical identification is correlated to 
the patterned response. The marketing term used often is 
"pattern recognition" . It can be a very powerful tool. 

[0061] In some cases an array of polymers that differ in 
physically properties, primarily solubility is used. The 
properties are measured electrically, change in resistance. The 
detector is about the size of a stamp. In other cases, an array 
of metal porphyrins is used. The pattern that develops upon 
exposure is read optically. Both of these cases require pattern 
recognition algorithms for chemical identification. 

Frequency Modulated Infrared Spectroscopy (FMIRS) 
[0062] Infrared spectroscopy is a group of techniques that 
measure the energy absorbed through molecular bond vibration and 
molecular rotation. Every chemical species has a unique IR 
fingerprint. The significance of these techniques is that 
chemical structural information is measured positively 
identifying the chemical when vibrational energies are collected 
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over sufficient range. IR and MS are the two techniques under 
this heading that provide positive specific chemical 
identification. The sensitivity is proportional to the number 
of molecules that interact with the energy source. The number 
of interactions is increased by increasing the distance the 
energy source travels through the sample. 

[0063] FMIRS is not a standard acronym in IR spectroscopy. 
There are a number of techniques in IR spectroscopy that use 
frequency or amplitude modulation. 

Biotechnology 

[0064] Microbiologists and geneticists have made great 
strides genetically engineering microbes that can specifically 
interact with chemical species, explosives, drugs, viruses, etc. 
The interaction results in some change in physical property of 
the supported microbe that can be measured. Most elegant of 
these is bioluminescence, biological light, emitted proportional 
to the interaction. Sufficient specificity is often attained 
for chemical characterization. A microbe must be developed for 
each chemical species of interest. 

Surface Acoustic Wave Resonator (SAW) 

[0065] A SAW detector is a crystal resonator held isothermal. 
A sample is passed over the detector. As chemical species in 
the air passing over the detector, they condense resulting in a 
change in resonator frequency. The amount of shift is a 
function of the number or concentration of chemical species in 
the air above the resonator. 
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SUMMARY OF THE INVENTION 



[0066] The present invention provides methods for 
interconnecting sub-systems, where each sub-system utilizes a 
single technology, and for identifying combinations of multiple 
technologies, where the combinations have a figure of merit 
satisfying a given criterion. 

[0067] The method of this invention for detecting the 
presence of a characteristic in an object includes: 

(a) providing three or more detection sub-systems, 

(b) first determining, utilizing one of the three or 
more detection sub-systems (detection technologies), 
whether the characteristic is present in the object 
whether the characteristic is present in the object, 

(c) subsequently determining, utilizing a subsequent one 
of the three or more detection sub-systems 
(detection technologies), whether the characteristic 
is present in the object, and, 

(d) finally, determining utilizing another one of the 
three or more detection sub-systems (detection 
technologies), if previous determinations for a 
majority of the three or more detection sub-systems 
are not identical, whether the characteristic is 
present in the object, and, 

(e) applying a predetermined criterion, after 
determinations by at least a majority of the three 
or more detection technologies, to detect the 
presence of the characteristic in the object. 

[0068] When it is desired to perform trade-offs and 
optimization of a system design (including systems capable of 
detecting the presence of more than one characteristic in an 
object), the method of this invention includes: 
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(a) provide a candidate system comprising at least two 
detection sub-systems; 

(b) performing a Monte Carlo analysis of detecting a 
presence of the at least one characteristic in an 
object; 

(c) calculating a pre-defined figure of merit for the 
candidate system; 

(d) selecting the candidate system if the figure of 
merit value satisfies a pre-determined criterion. 

[0069] Systems obtained by the above methods are also 
disclosed. 

[0070] For a better understanding of the present invention, 
together with other and further objects thereof, reference is 
made to the accompanying drawings and detailed description, and 
its scope will be pointed out in the appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0071] Figure 1 depicts a flowchart of an embodiment of a 
method of this invention; 

[0072] Figure 2 depicts a Block/Flow diagram for a subsystem 
for determining the presence of a characteristic in an object; 
[0073] Figure 3 illustrates the Probabilities set for the 
subsystem of Figure 2 ; 

[0074] Figure 4 depicts a Block/Flow diagram for a Level 1 
system (Assuming all levels have PD= 90% and PFA = 10%); 
[0075] Figure 5 depicts a Block/Flow diagram for a Level 2 
system (Assuming all levels have PD= 90% and PFA = 10%); 
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[0076] Figure 6 depicts a Block/Flow diagram for an 
embodiment of a system of this invention (Assuming all levels 
have PD= 90% and PFA = 10%); and, 

[0077] Figure 7 depicts a flowchart of another embodiment of 
the method of this invention. 

DETAILED DESCRIPTION OF THE EMBODIMENTS OF THE INVENTION 

[0078] The methods and systems of the present invention 
implement interconnected systems for detecting a presence of a 
characteristic in an object. It should be noted that "object" as 
used herein refers to the situation in which the characteristic 
is being detected. An "object" as used herein below may be, but 
is not limited to, an object, a patient, a region or location on 
a body, an object at a time and place, or the interior of an 
object. A "characteristic" as used herein refers to that which 
is being detected. A "characteristic" as used herein below may 
be, but is not limited to, an undesired substance (such as, but 
is not limited to, an explosive) within an object, whether an 
undesired substance is not present in an object, a particular 
disease at a location in a patient, location of targets, 
obstacles and objects in an area or volume. 

[007 9] Strategies for interconnecting sub-systems, each sub- 
system utilizing a single detection method, and for selecting 
combinations of sub-systems designed according to one of the 
simpler strategies are disclosed hereinbelow. 

[0080] Figure 1 depicts a flowchart of an embodiment of a 
method of this invention. Referring to Fig. 1, three or more 
detection (N) sub-systems are provided (step 20, Fig. 1). The 
first sub-system is utilized to determine whether the 
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characteristic is present in the object or the characteristic is 
absent in the object (step 30, Fig. 1). Subsequent sub- systems, 
up to the last sub-system, are utilized to determine whether the 
characteristic is present in the object (steps 30, 40, Fig. 1). 
After a majority of objects (>N/2) of the sub-systems have been 
utilized to determine whether the characteristic is present, if 
previous determinations for a majority of the three or more 
detection sub-systems are not identical, another subsystem is 
utilized to determine whether the characteristic is present in 
the object(step 50, Fig. 1). The process continues until 
determinations for a majority (>N/2) of the three or more 
detection sub-systems are identical or until all sub-systems 
have been utilized (step 55, Fig. 1). If identical 
determinations for a majority (>N/2) of the three or more 
detection sub-systems are obtained, a final determination of the 
presence of the characteristic is obtained (step 60, Fig. 1). If 
a majority of identical determinations are not obtained and all 
sub-systems have been utilized (as, for example, in an 
embodiment utilizing an even number of sub-systems), the final 
determination of the presence of the characteristic is obtained 
utilizing a pre-determined criterion. 

[0081] A detection technology can be represented by the flow 
diagram in Figure 2. Referring to Fig. 2, when an object is 
examined, subsystem 1 produces an Alarm or OK response if the 
object does or does not contain a given characteristic, 
respectively. The probabilities of correct responses and errors 
are shown in Figure 3 . 

[0082] Referring to Fig. 3, PD, the probability of detection, 
is the probability of detecting a characteristic if there is a 
characteristic present. PFA, the probability of false alarm, is 
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the probability of detecting a characteristic , if there is none. 
This leads to threat mitigation procedures. 

[0083] Referring again to Fig. 3, 1-PD is the probability of 
failing to detect a characteristic. In the case where the 
characteristic is an explosive, a failure to detect, when there 
is an explosive present, could lead to an explosion. 

[0084] If PFA is high it leads to costly time consuming 
threat mitigation procedures. If PD is low and 1-PD is high 
this could lead to failure to detect the characteristic. In the 
case where the characteristic is an explosive, a failure to 
detect could lead to an explosion, loss of life, and extreme 
costs. Note that PD and PFA do not have to add to 1. Tuning 
any single system (or technology) tends to raise or lower both 
PD and PFA, together. 

[0085] Combinations of technologies can be considered to 
raise the system probability of detection (SPD) and lower the 
system probability of false alarm (SPFA), but not all strategies 
produce improvements. Multi-technology systems require 
independent technologies in the subsystems. 

[0086] The following example strategies for a Level 1 System, 
a Level 2 Alarm Mitigation System, three level decision system, 
and the majority system assume PD= 90% and PFA = 10% for all 
technologies. Figures 4-6 show the result of examining 100 
objects known not to have the characteristic (without 
explosives, for example) and 100 objects containing the 
characteristic (with explosives, for example). The System 
Probability of Detection (SPD) and the System Probability of 
false alarm (SPFA) can be calculated. 
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[0087] The Level 1 System 100, shown in Figure 4, has a 
baseline SPD = 90% and a baseline SPFA = 10%. 

[0088] Figure 5 shows a 2 level Alarm Mitigation System 
strategy. The combined system improved the SPFA to 1% at the 
expense of significantly reducing SPD to 81%. 

[0089] Figure 6 shows a 3 level Majority System strategy. 
The system 200 is structured to require either two or three 
technologies to determine that the characteristic is not present 
("OK the object") or two of three technologies to Alarm the 
object before a determination is made. The strategy is 
effectively having each the technologies vote, and the third 
level technology is used as the tiebreaker. In practice, the 
level 2 and level 3 technologies would probably be slower than 
the level 1 technology. The system could actually derive a 
significant speed-up using a state variable scheduler to direct 
objects to the least busy level. This is the only strategy that 
can improve both the SPD and SPFA at the same time and can 
improve the throughput. 

[0090] In one embodiment, where only three detection 
techniques (sub-systems) are utilized, the probability of 
detection (P D1 ) of the sub-system that has the highest 
probability of detection (110, Fig. 6) satisfies the following 
relation: 

(1- P D1 ) * [(1- P D i)(l- P d2 ) + (1- P D1 )P D2 (1- P D3 ) 

+ P D i(l- Pd2)(1" P D3)1/ 
where P D2 is a probability of detection of a second detection 

sub-system (120, Fig. 6), 
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P D3 is a probability of detection of a third detection sub-system 
(130, Fig. 6). 

[0091] In another embodiment, the probability of false alarm 
(P FA1 ) of the sub-system (also known as the probability of false 
positives) that has the lowest probability of false alarm (in 
one embodiment, 110, Fig. 6) satisfies the following relation: 

Pfai * [ (Pfai)(Pfa 2 ) + (1 - Pfai)(Pfa 2 )(Pfa3) 
+ (Pfai)(1 - P F a 2 )(Pfa3)] where, 

P FA2 is a probability of false alarm of a second detection sub- 
system (120, Fig. 6), and 

P FA3 is a probability of false alarm of a third detection sub- 
system (130, Fig. 6). 

[0092] Embodiments in which both of the above conditions are 
satisfied by the same sub-system are also possible. 

[0093] In the above embodiments, the probability of detection 
for the system and the combined probability of false alarm for 
the system are jointly improved over that of a single detection 
method (technology). It should be noted that the above described 
conditions can be met if all sub-systems (technologies) have the 
same probability of detection and the same probability of false 
alarm and the probability of detection is greater than 0.5 and 
the probability of false alarm is less than 0.5. It should also 
be noted that the sub-system that has the highest probability of 
detection is not, in every embodiment, the same sub-system that 
has the lowest probability of false alarm. It should be further 
noted that although, in the embodiment in which the same sub- 
system has the highest probability of detection and the lowest 
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probability of false alarm shown in Fig. 6, although that sub- 
system is shown as the first sub-system in Fig. 6, the sub- 
system that has the highest probability of detection and/or the 
sub-system that has the lowest probability of false alarm could 
be placed in a different order (it could be, but not limited to, 
the second sub-system) . 

[0094] In another embodiment, the number of detection sub- 
systems is an even number. In that embodiment, a predetermined 
criterion is utilized to detect the presence of the 
characteristic. If a majority (>N/2) of identical determinations 
is obtained, the presence of the characteristic is detected. One 
embodiment of the predetermined criterion is, if a majority 
(>N/2) of identical determinations is not obtained, detecting 
the presence of the characteristic if half of the detection sub- 
systems determine that the characteristic is present. 

[0095] In the embodiment in which the detection of the 
presence of a characteristic in an object is the detection of 
explosives inside objects, each one of the detection sub-systems 
(110, 120, 130, Fig. 6) can utilize, but is not limited to, one 
of the following technologies: X-ray technologies, Neutron-Based 
Bulk Explosives Detection Technologies, Nuclear non-neutron 
based Bulk Explosive Technologies, Nuclear Quadrupole Resonance 
Technologies, and Trace Detection Technologies. The Nuclear non- 
neutron based Bulk Explosive Technologies include Gamma 
Resonance Technology (GRT) and Gamma Radiography. It should be 
noted that other technologies (such as future technologies) can 
be utilized if available. 

[0096] The Majority Strategy described above (and in Fig. 6) 
can be used to predict the functionality of technology 
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combinations. The probability of detection and probability of 
false alarm represent a base line for combined technology 
systems. However , for a large number of sub-systems or for 
detecting presence of more than one characteristic in an object, 
other methods can be used. The Monte Carlo method can be used to 
determine the probability of detection or false alarm for any 
combination of sub-systems and for detecting presence of more 
than one characteristic in an object. The Monte Carlo method is 
a computational technique that uses random numbers to simulate 
the behavior of a system. By conducting a large number of 
simulation objects, properties of the system, which would not be 
obtainable by derivation from closed form mathematical solution, 
can be ascertained. An example, but not limited to, of detecting 
the presence of more than one characteristic in an object is the 
detecting of several particular explosives. 

[0097] A flowchart of an embodiment of the method of this 
invention utilizing the Monte Carlo method is shown in Figure 7. 
Referring to Fig. 7, a candidate combination including two or 
more detection sub-systems is selected (step 210, Fig. 7). A 
Monte Carlo analysis of detecting the presence of the one or 
more characteristics in the object for the system including the 
two or more detection sub-systems is performed (step 220, Fig. 
7). Utilizing the results of the Monte Carlo analysis and system 
properties, a pre-defined figure of merit (also referred to 
herein below as a "cost") is calculated (step 230, Fig. 7). If 
other sub-system combinations are available for analysis, the 
previous steps are repeated (step 240, Fig. 7). Those 
combinations having a Figure of Merit satisfying a pre- 
determined criterion (in one embodiment, larger than a given 
value) are selected (step 250, Fig. 7). 
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[0098] In one embodiment, three sub-systems are selected 
utilizing the 3 level Majority System strategy of Figure 6. In a 
particular embodiment of the Majority System strategy, the 
probability of detection (P D1 ) of the sub-system that has the 
highest probability of detection (110, Fig. 6) satisfies the 
following relation: 

(i- p D1 ) * [(i- p D1 )(i- p D2 ) + (i- p D1 )p D2 (i- p D3 ) 
+ p D1 u- p D2 )(i- p D3 )] 

where , 

P D2 is a probability of detection of a second detection sub- 
system (120, Fig. 6), 

P D3 is a probability of detection of a third detection sub-system 
(130, Fig. 6). 

[0099] In another embodiment, the probability of false alarm 
(P FA1 ) of the sub-system that has the lowest probability of false 
alarm (in one embodiment, 110, Fig. 6) satisfies the following 
relation: 

Pfai * [ (Pfai)(Pfa 2 ) + (1 - Pfai)(Pfa 2 )(Pfa3) 
+ (P FA1 )(1 - P F a 2 )(Pfa3)] where, 

P FA2 is a probability of false alarm of a second detection sub- 
system (120, Fig. 6), and 

p fa3 is a probability of false alarm of a third detection sub- 
system (130, Fig. 6). 

[00100] Embodiments in which both of the above conditions are 
satisfied by the same sub-system are also possible. 

[00101] An embodiment of the method of this invention, 
utilizing the Monte Carlo method, for performing design trade- 
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offs, selecting candidate sub-systems according to a pre-defined 
cost function (Figure of Merit) where the detection of the 
presence of a characteristic in an object is the detection of 
explosives inside objects is described herein below. 

[00102] Referring again to Fig. 5, assume further that a 

probability of an explosive being present (PE) in 1 out of 
100,000,000 objects or 10" 8 and a probability of non-explosives 
being in the balance of the objects (PNE) of 0.99999999% 

[00103] The cost component (CC) consists of the sum of the 
cost of failure to detect an explosive and the cost of false 
alarms. Multiplying each of these by the cost of the response 
and the number of objects per year gives the expected system 
costs component per year. Equation 1 gives the cost component 
per year. 

Eq. 1. CC = (N*PE*(1-SPD)*CE) + (N*PNE*SPFA*CFA) 

where CC is the system cost component, $/year; N is the number 
of objects per year being screened; SPD is the system 
probability of detection; CE is the cost of an explosive not 
being detected, $/object; SPFA is the system probability of 
false alarm; and CFA is the cost of a bomb squad examining a 
false alarm, $/object. 

[00104] If it is assumed that the cost of an explosion is 
$1,000,000,000, the cost for a bomb squad handling a false alarm 
is $1,000, and the number of screenings is 100,000,000 objects 
per year, then expected cost component, CC, is 

Eq. 2. CC = $1,200,000,000 / year 
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[00105] It may be possible to optimize the component cost, CC 
by examining the derivative of CC with respect to SPD and SPFA. 

[00106] Other cost components include, but are not limited to, 
purchase price, installation, maintenance, operating cost. 

[00107] In the embodiment of the application of the method of 
this invention to selecting candidate combinations of sub- 
systems for an Explosive Detection System (EDS) presented below, 
scoring criteria was utilized before applying the method of this 
invention and a ranking of individual EDS sub-systems was 
performed utilizing the scoring criteria prior to applying the 
method of this invention. Although the scoring criteria and the 
preliminary ranking are not necessary to perform the method of 
this invention, the scoring criteria and the preliminary ranking 
are described below since reference is made to them in the 
embodiment of the application of the method of this invention 
disclosed herein below. 

EDS Scoring Criteria 

[00108] The Evaluation Criteria used in the study are: 
Probability of Detection, Specificity, Probability of False 
Alarm, Unit Cost, Annual Maintenance Cost, Staff Requirements, 
Availability to Market, Footprint area, Weight, Throughput, 
Automation level and Defeat Immunity. 

[00109] The Scoring Criteria, which establish scores based on 
evaluation criteria performance and corresponding Weighting 
Factors, which establish the relative importance of Evaluation 
Criteria are shown in Table 3 . 
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Table 3: EDS Scoring Criteria of this invention 



EDS Scoring Criteria 



Rank 


Best 






Worst 














Weight 
Factor 

note 10 


Score 


10 


9 


8 


7 


6 


5 


4 


3 


2 


1 


Evaluation Criteria 






















Prob. Of Detection, % 


>99 


99 - 98 


98 - 97 


97 - 95 


I 95-90 


90-85 


85-80 


80-75 


75-70 


<70 


20 


Specificity (number of material 
parameters measured) 


5 or more 


4 


3 


2 














14 


Prob. Of False Positive, % 


<0.1 
< 200 


0.1 -0.2 0.2-0.3 0.3-0.5 


0.5-1 


1 -5 


5-10 j 


10-20 


20 - 30 


> 30 


14 


Unit Cost . $K 


200 - 500 


500 - 800 


800-1000 


1000- 1200 


1200-1500 


1500 - 2000 


2000 - 3000 


3000 - 5000 


>5000 


13 


Annual Maintenance Cost, $K 


< 20 


20-50 


50-80 


80-100 


100-120 


120-150 


150-200 


200 - 300 


300 - 500 


>500 


2 


Annual Operational Support , Full 
Time Equivalent Staff 


0 


0.5 


1 


1.5 


2 


2.5 


3 


4 


5 


>5 


2 


Availability to Market, months 


Now 


1 -6 


6-12 


12-24 


24-36 
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Footprint, Sq.-Ft. 


<50 


50-100 


100-200 


200 - 500 


500- 1000 


1000- 1500 


1500-2000 


2000 - 2500 


2500 - 3000 


>3000 


6 


Total Weight, lb 


< 500 


500 - 1000 


1000- 1500 1500 - 2000 


2000 - 5000 


5000- 10000 


1 0000- 1500< 


5000 - 2000 


20000 - 25000 


>25000 


3 


Throughput (objects/mi n) 


> 25 


25-15 


15-5 


5-1 














3 


Automation 


Auto level 1 
and 2 


Auto level 1 
and manual 
level 2 


Auto level 1 


Manual level 
1 














1 


Defeat Immunity 


Can't be 
defeated 


Defeatable 
by expert , 
indicates 
defeated 


Defeatable 
easily and 
inciateds 
defeated 


Defeatable 
and does not 
indicate 
defeated. 
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EDS Scores— 

[00110] Scores can be obtained for various explosive detection 
sub-systems (EDS). The characteristics of each EDS, recorded on 
an EDS Analysis Factors Form, are the basis for the scores. 

[00111] The Figure of Merit is the sum of the products of the 
scores multiplied by the corresponding weighting factor, all 
divided by 100. An example of an EDS Ranking is shown below in 
Table 4 . 
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Table 4: EDS Ranking of this invention 



EDS Scores 



Technology 


DE 
CT 


Fixed 
Emitter/ 
Detector 


Diff 
Xray 
1 


Diff 
Xray 
2 


Diff 
Xray 
3 


Side 
Scan 
Xray 


Fixed 
Emitter/ 
Detector 


LPFNA 
1 


SPFNA 


LPFNA 
2 


TNA 


NQR 
1 


NQR 

2 


Prob. of Detection 


3 


2 


4 


3 


4 


2 


2 


4 


4 


3 


3 


3 


3 


Specificity, note 5 


3 


3 


4 


4 


4 


3 


3 


3 


4 


3 


1 


4 


4 


Prob. of False Positive 


3 


2 


4 


4 


4 


2 


2 


4 


4 


4 


3 


4 


4 


Unit Cost 


2 


2 


2 


2 


1 


3 


2 


3 


1 


2 


3 


4 


4 


Annual Maint. Cost 


4 


1 


1 


1 


1 


2 


2 


3 


1 


3 


4 


4 


4 


Annual Operational Support 


3 


3 


4 


4 


4 


3 


3 


4 


4 


4 


4 


4 


4 


Availability to Market 


2 


3 


3 


2 


2 


1 


1 


2 


4 


1 


3 


2 


2 


Footprint 


3 


3 


3 


3 


2 


4 


3 


2 


1 


2 


3 


4 


4 


Total Weight 


1 


1 


1 


1 


1 


3 


1 


1 


1 


2 


2 


3 


3 


Throughput, objects/minute 


3 


4 


1 


2 


1 


2 


4 


2 


3 


1 


1 


2 


1 


Automation Level 


3 


3 


2 


2 


2 


3 


3 


2 


4 


2 


2 


2 


2 


Defeat Immunity 


3 


3 


1 


1 


1 


3 


3 


4 


4 


4 


4 


1 


1 


Figure of Merit (note 11)= 


2.67 


2.4 


2.94 


2.67 


2.65 


2.41 


2.22 


3.14 


3.17 


2.71 


2.7 


3.11 


3.08 


notes 






.4 


4,8 


4,8,12 


8,9 












4 


4,6,7 



NOTES : 

1. Objects = tubs with flats 18 R x13"x1 1"H, sacks with objects 40"x30" flat or 36"Lx24"D, 
Outiiers<107"L. 

2. Total Weight includes device + shielding (if applicable) 

3. Floor space is the foot print area including shielding, square feet 

4. Detect crystals only, therefore may not detect dynamite. 

5. Specificity scoring: whole object at once=0pt; 2D=1pt; 3D=1pt; density=1pt; eff Z=1pt; 
linear x-ray attenuation coefficients pt; measure N=2pt; measure O, H, C=2pt; molecular 
ID=1pt; all scores are additive. 

6. Based on TSA and FAA test estimates. 

7. Aluminum foil will defeat NQR. 

8. Ge detector requires LN or cooling. Can see crystals only so can't see Dynamite. 

9. The Weight Factors total 100%. 

[00112] The EDS Ranking is useful for comparing the individual 
explosive detection systems over a wide range of evaluation 
criteria. 

[00113] A Composite Systems Comparison was performed utilizing 
the method of this invention for selecting candidate 
combinations of sub-systems as shown in Fig. 7. The Composite 
Systems Comparison identifies which of the EDS sub-systems 
produce combined systems in which the Figure of Merit indicates 



10. The Figure of Merit is the sum of each score multiplied by the weighting 
factor, all divided by 100. The Best is 4 and the worst is 1 . 

11. This price based on unit cost. No quantity purchase data available. 

12. Operational Support is on a basis of one shift per day. 
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that the performance is above a given performance value and the 
"cost" is below a given "cost" value (with a goal of the best 
performance for the lowest cost). The input data are the values 
from the Analysis Sheets or the midrange of the EDS Ranking. 
Technologies that cannot detect noncrystalline (NC) materials 
are assumed to have a PD NC of 30%. 

[00114] The EDS sub-systems from the EDS Ranking study are 
combined in groups of one, two, and three, using the Majority 
Strategy described herein above. Each combination is analyzed 
by using a Monte Carlo simulation to determine its performance. 
The Monte Carlo analysis interprets a random numbers generator 
to determine if the object being screened is: 

Does not contain an explosive 

Contains a crystalline explosive 

Contains a noncrystalline explosive 

[00115] Some of the technologies do not detect noncrystalline 
explosives well. Thus, separate probabilities of detection are 
used as inputs to the calculations. Actual data on how much 
crystalline and noncrystalline explosive are present in objects 
is not known. The Bureau of Alcohol, Firearms, and Tobacco 
(ATF) website reports that 86% of the explosives stolen each 
year are crystalline explosives (dynamite). Detecting an 
explosive as common as dynamite should be an important 
consideration . 

[00116] The Monte Carlo analysis outputs are: 

(a) PD NC Probability of detection for noncrystalline 

materials. The value is from the Analysis Sheet or the 
midrange of the EDS Ranking. The value used was assumed to 
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be 30% for technologies that cannot detect noncrystalline 
materials 

(b) PD C Probability of detection for crystalline materials. 

(c) SPD Probability of detection for the combined system. 

(d) Specificity Combines the specificity of all input 
technologies . 

(e) SPFA Probability of false alarm for the combined system. 

[00117] Unit Cost Adds the cost of all combined sub-system. 

(a) Maintenance Cost adds the maintenance cost of all combined 
sub-systems . 

(b) Staff Adds the number of full time people required for 
each subsystem per shift. 

(c) Availability Months until availability to the market. 

(d) Floor Space Adds the area required for each subsystem. 

(e) Weight Adds the weight of each subsystem. 

(f) Throughput Adds the time required for each subsystem to 
screen the object. 

(g) Automation Combines the automation levels for each 
subsystem. 

CFOM Cost Figure of Merit. (See below.) 

[00118] The CFOM was derived by calculating the actual cost 
per year for each combination of sub-systems , using the 
equation: 

CFOM = 120*(1-PD) + 120*PFA + unit cost + maint. Cost + 
staff cost (all in $ millions) 

[00119] This equation was derived by estimating that the cost 
for an explosive detection error to be $120M/year times the 
probability of a missed detection and the cost of false alarms 
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is $120M/year times the probability of false alarm. The cost of 
an explosion and the cost of false alarms are kept similar so 
the Monte Carlo analysis wouldn't overly emphasize one or the 
other. 

[00120] There are a number of methods of selecting EDS sub- 
systems as candidate sub-systems for the design of a system. It 
should be noted that sometimes different methods lead to the 
same conclusions. 

[00121] One method of analysis could be, but is not limited 
to, selecting an acceptable SPD and SPFA and then eliminating 
the EDS sub-systems with the lowest score. 

[00122] Another method of analysis could be, but is not 
limited to, selecting an acceptable SPD and SPFA and then sort 
by CFOM and eliminate the most expensive combinations. 

[00123] Alternately, the analysis could be done, but is not 
limited to, by selecting an acceptable SPD and SPFA and then 
sort by CFOM as before, and the sub-systems that appear most 
frequently in the most expensive combinations could be 
eliminated. For example, for SPD>90% and SPFA<1%, sub-systems 
that appear most frequently in the most expensive combinations 
could be eliminated from being candidate sub-systems. 

[00124] Applying the analysis methods described above, sub- 
systems can be selected for a desired outcome. 

[00125] It should be noted that, although EDS scores were 
presented for a group technologies, other technologies can be 
utilized. The system of this invention can utilize, but is not 
limited to, one of the following technologies: X-ray 
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technologies, Neutron-Based Bulk Explosives Detection 
Technologies, Nuclear non-neutron based Bulk Explosive 
Technologies, Nuclear Quadrupole Resonance Technologies, and 
Trace Detection Technologies. The Nuclear non-neutron based Bulk 
Explosive Technologies include Gamma Resonance Technology (GRT) 
and Gamma Radiography. Other technologies (such as future 
technologies) can be utilized if available. 

[00126] It should also be noted that although in the 
embodiments presented above the detection of the presence of a 
characteristic in an object is the detection of explosives 
inside objects, other embodiments of characteristics present in 
objects are also possible and within the scope of this 
invention. 

[00127] Although the invention has been described with respect 
to various embodiments, it should be realized that this 
invention is also capable of a wide variety of further and other 
embodiments all within the spirit and scope of the appended 
claims • 

[00128] What is claimed is: 
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